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Abstract We present the first high-spatial resolution 
near-infrared (NIR) imaging of NGC 604, obtained with 
the NICMOS camera aboard the Hubble Space Tele- 
scope (HST). These NICMOS broadband images reveal 
new NIR point sources, clusters, and diffuse structures. 
We found an excellent spatial correlation between the 
8.4 GHz radio continuum and the 2.2 /xm nebular emis- 
sion. Moreover, massive young stellar object candidates 
appear aligned with these radio peaks, reinforcing the 
idea that those areas are star-forming regions. Three 
different scaled OB associations are recognized in the 



NICMOS images. The brightest NIR sources in our 
images have properties that suggest that they are red 
supergiant stars, of which one of them was previously 
known. This preliminary analysis of the NICMOS im- 
ages shows the complexity of the stellar content of the 
NGC 604 nebula. 
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1 Motivation for a detailed study of NGC 604 

Giant H II regions (GHRs) are among the most lumi- 
nous objects that can be individually identified very 
distant galaxies. In these regions star formation oc- 
curs at extremely high rates, hence they are also known 
as starburst regions. Given the young age and large 
number of massive stars in GHRs, some of them are 
still found embedded in their dusty parental molecular 
clouds. The morphology of these regions changes quite 
rapidly during the first few million years after the first 
massive generation is born. Excellent pictures showing 
these evolutionary morphologies are the H II regions in 
the irregular galaxy N GC 4214 (jMaiz Apellaniz! bOQOl : 
MacKentv et alj|2000( ). 

The natural extension to larger scales are the star- 
burst galaxies which can be observed at cosmological 
distances and, therefore, can be used as astrophysical 
signposts to trace the star formation history of the Uni- 
verse. However, due to their large distances, we must 
rely only on their major global properties, such as the 
behavior of the strong recombination lines originating 
in the ionized nebula. In order to improve our knowl- 
edge of these objects in distant galaxies, it is necessary 
to study resolved nearby examples. 30 Doradus in the 
LMC and NGC 604 in M33 are the two largest GHRs in 
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Fig. 1 2MASS Tfs-band image of the region of NGC 604. 
IR sources labeled as 1, 2, 3, and 4 are associated directly 
with the molecular clouds (contours, adapted from Engargi- 
ola et al. 2003) . Fields observed with the NICMOS camera 
2 are marked as rectangles in full lines. North is up and 
East to the left. The point source labeled as "F" is a galac- 
tic foreground star. 



the Local Group (LG), so they are the prime candidates 
for such study. 



NGC 604 is located at a distance of 840 kpc (jFreedman 
2001). At such a distance, 1" is equivalent to 4pc 



into three sets: (a) slitless objective-prism FUV spec- 
troscopy obtained using ACS/SBC, (b) multi-filter (6 
broadband and 2 narrowband) NUV to optical imag- 
ing obtained using ACS/HRC, and (c) NICMOS/NIC2 
broadband imaging. These datasets were obtained un- 
der HST proposals 10419 and 10722. The existing data 
include STIS-NUV objective-prism spectroscopy and 
photometry obtained under HST/GO program 9096, as 
well as archival HST/WFPC2 imaging. Additional ob- 
servations were obtained using Gemini North/NIRI in 3 
broadband and 3 narrowband NIR filters and with the 
Combined Array for Research in Millimeter-wave As- 
tronomy (CARMA) in multi-configurations C, D, and 
E, to observe the CO(l — > 0) transition with a velocity 
resolution of 3kms _1 , and a beam size ranging from 
3"-5". 

The main goals of this project are: to obtain spec- 
tral classifications for about 200 stars, measuring their 
spectral energy distribution (SED) from 130 nm to 
2.2 /im, to identify embedded very young stellar pop- 
ulations hidden inside dust clouds, to measure the 
extinction law and its possible variations as a func- 
tion of the environmental conditions, and to analyze 
the relationship between the hot stars and the sur- 
rounding gas. The combination of all data will al- 
low to obtain the most complete and detailed extinc- 
tion corrected Hcrtzsprung-Russell diagram of a re- 



tion corrected rlcrtzsprung-ttusscil diagram ot a re- 
et j&y id SOB A. T he SEDs obtained w ill be analyzed 



and consequently the capabilities of HST are essen- 
tial to study the region with high spatial resolution. 
NGC 604 is powered by a massive young cluster without 
a central core that con tains over 200 O and WR stars 
( Maiz Apellania l200ll and references therein). This 



GHR is the best nearby exampl e of a scaled OB associ- 
ation, or SOBA, as defined bv iMaiz Apellanizl (|200lh . 
a more extended type of object than 30 Dor, which is a 
Super Star Cluster (SSC). The combination of similar- 
ities and differences between those objects makes the 
combined analysis of both 30 Doradus and NGC 604 
a necessary step for the creation of a template for the 
understanding of far-away starbursts. 

In this contribution we describe the gathered obser- 
vations and general objectives, and we present a pre- 
liminary analysis of the new HST/NICMOS images. 



2 Observations 

We are dedicating a large observational effort to do a 
thorough multi-wavelength study of NGC 604, from the 
far ultraviolet to the near infrared with additional high- 
resolution millimetric observations, combining existing 
data with the new ones. The new HST data are divided 



with CHORIZOS |Maiz Apellanizll2004h as the primary 
code. 



3 A NICMOS view of NGC 604 

Six HST orbits were requested and allocated to ob- 
serve six selected fields in and near NGC 604 with the 
NICMOS camera 2 (NIC2), using the broadband filters 
F110W, F160W and F205W (similar to the J, H, and K 
passbands, respectively). Figure [1] shows the observed 
fields on a Ks image extracted from 2MASS. A 37" x 32" 
(~ 140 pc) area centered on the core of NGC 604 was 
observed in four tiles using a dithering strategy. An ad- 
jacent field to the south of the core, centered on a bright 
NIR source was also observed. In addition, sky frames 
were obtained in a nearby (1' = 240 pc away) field for 
background subtraction and to take into account the 
stellar population of the M33 disk in the subsequent 
photometric analysis. 

Molecular clouds associated with NGC 604 have been 
detected in CO by lEngargiola et al. (|2003l and refer- 
ences therein). The 2MASS ifs-band image (Fig.[T|) 
shows a concentration of point-like sources in the core 
of NGC 604, some of them surrounded by an extended 



NICMOS view of NGC 604 
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Fig. 2 NICMOS F110W (top panel) and F205W (middle panel) mosaics of NGC 604. The bottom panel is an explanatory 
diagram for some structures discussed in the text. SOBAs labeled as "1", "2" and "3" are the main SOBA, the secondary 
SOBA and a third one located to the south, respectively. The "F" letter indicates the optically visible gaseous filame nts. 
Flaps are areas with high optical extinction. "R" refers to the red supergiant star discovered bv lTerlevich et all l)l99rt l 
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halo. In contrast, the "background" field seems clean of 
such bright sources. The 2MASS IR sources are associ- 
ated directly with the molecular clouds. This morpho- 
logical relationship between infrared sources and molec- 
ular clouds resembles the 30 Dor Nebula. 

The new NICMOS broadband imaging of NGC 604 
reveals new IR point sources, clusters and diffuse struc- 
tures, improving the spatial resolution of the 2MASS 
data by more than an order of magnitude. Figure [2] 
shows the F110W and F205W mosaics of five NIC2 
fields in the observed area of NGC 604. The variety 
of observed structures is striking. The two optically 
detected SOBAs (main and secondary) are easy to dis- 
cer n. They are encircled by la beled ellipses in Fig. [2] 
(see Maiz Apellaniz et al.lbooi] for detailed explanation 
of the optical structures). An additional SOB A af- 
fected by larger extinction is present to the south (la- 
beled as "SOB A 3" ). The two cavities associated to the 
SOBAs are defined. Cavity A (associated to the main 
SOBA) is limited by a bright rim to the south. This 
rim is brighter in the F205W filter, when compared 
with the F110W filter, indicating that it is affected by 
a large extinction. The extended nebular emission to 
the southwest in the F205W im age is coincident with 



the structure labeled as "flap" bv lMaiz Apellaniz et al 



(2004J), an (optical) obscuration area which is a dusty 
screen associated with the molecular cloud in that re- 
gion. This main CO clo ud lies on the southern ed ge of 
the bright nebular rim. IChurchwell &: Gosd (] 19991 ) and 
Maiz Apellaniz et al. (|2004f ) have analyzed the extinc- 
tion in NGC 604 using radio continuum, Ha, and H/3 
images. They have found a very good correlation be- 
tween the CO emission and the total optical extinction, 
which indicates that the dust is associated with molec- 
ular gas. Some of the H II gas is located behind large- 
optical-depth well-defined dust "flaps" that make it in- 
visible in the optical. Figure [3] s hows the 8.4 GHz radi o 
continuum contours (from (jChurchwell fc Gosslll999l n 
superimposed on the F205W image. Again, the cor- 
relation between radio and 2.2/im emission seems to 
be very good. The radio peaks C and D are super- 
posed to the southwest NIR nebular emission. The ra- 
dio peaks A and B are coincident with nebular knots. 
Peak A is one of the bri ghtest Ha structures (see 
( Maiz Apellaniz et al.l [2004] ) ). and the same happens 
in the NIR. Gemini/NIRI images (not shown) demon- 
strate that most of such emission is produced in the 
warm gas as Brackett 7 emission line. A very interest- 
ing result is the remarkable spatial correlation between 
the massive young stellar object candidates and the ra- 
dio peaks, which suggests that those areas are prime 
candidates for massive star forming regions. 




Fig. 3 Contour diagram of the 8.4 GHz radio continuum 
(adapted from Churchwell & Goss 1999) superposed on the 
F205W NIC2 image. Labels indicate the radio continuum 
structures (see Mai'z Apellaniz et al. 2004 for an explana- 
tion). Circles mark the position of the massive young stellar 
object candidates 



4 NICMOS photometry 

Point spread function photometry was performed in the 
NIC2 mosaics and "background" field. In total, we have 
more than 600 detections in the three filters with pho- 
tometric errors less than 0.3 mag. The photometric 
sensitivity is about 21.9 mag in the F205W filter, cor- 
responding to an unreddened early-B type star at the 
distance of M33. FigurcH] shows the color-color and 
color-magnitude diagrams for all sources brighter than 
^F205W = 19.5 (03 V star at M33). We decided to 
preserve our data in the NICMOS photometric sys- 
tem in order to avoid additional problems with color 
transformations to ground-based photometric systems. 
Main-sequence and red giant stars loci are plotted as 
continuum lines in the color-color diagram in Fig.fJJ 
They are synthesized using CHORIZOS, Kurucz mod- 
els, and NIC2 filter transmission curves. Reddening 
vectors (Ay = 16) for an early-type star and a red giant 
have been derived in the same way, adopting a normal 
reddening law. In the color-color diagram is possible 
to discern three main population of objects: (a) blue, 
hot stars with small o moderate reddening close to the 
main-sequence locus and along the reddening track; (b) 
red stars in the red giant locus; (c) a population of stars 
with infrared excess to the right of the early-type star 
reddening vector. 

The distribution of sources in the color-magnitude 
diagram shows that the blue sources (group a) cor- 
respond mostly to luminous hot stars associated to 
SOBAs (massive stars). Meanwhile, stars in the group 
(c) are mostly massive young stellar object (MYSO) 
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Fig. 4 NIC2 color-color (left) and color-magnitude (right) diagrams. This plot includes all sources with mF205w < 19.5. 
Color-color diagram: the main-sequence locus and the cool-giant branch are indicated as solid lines; reddening tracks for 
Av = 16 are also plotted. Color-magnitude diagram: red supergiant star locus is indicated as a dotted line in the upper 
part of the diagram; the asymptotic giant branch and red giant branch are also plotted as full lines. Symbols: "stars" are 
MYSO candidates; "squares" are hot stars; "dots" are mostly red giant stars; "black circles" are red supergiants stars 



candidates and/or extreme emission line objects. Ob- 
jects in group (b) are stars belonging to the red gi- 
ant branch, and they are more or less smoothly dis- 
tributed in the mosaic. The brightest objects in the 
field are five stars in group (b). Their positions in the 
color-magnitude diagram suggest that they are red su- 
pergiant stars. In fact, one of them has been spectro - 
scopically identified as such by iTerlevich et al.l (|1996[ ). 
This star is a close companion to a compact mas- 
sive clust er and two WR sta r ca ndidates (WR11 and 



WR7, see lDrissen et al.lll993l and lMaiz Apellaniz et al 
2004). The other four red supergiant star candidates 



are also located very close (few arc-seconds) to bona- 
fide WR stars, WR candidates or massive hot stars. 
The spatial neighborhood observed between the RSG 
stars and the early-type objects addresses the question 
of whether they belong to the same star generation or 
they just represent a superposition of different popula- 
tions, as it is the case of 30 Dor. In the latter GHR, 
Walborn fc Blades! (|1997I ) proposed the coeval existence 
of five distinct populations, nicely distinguished by age 
and distribution. Likewise, this issue indicates the com- 
plexity of the stellar content in a GHR. 

Among the MYSO candidates, a particular source 
deserves special consideration, owing to its large NIR 
colors. This object, labeled as "Y" in Figfig:radio and 
Fig.His barely detectable in the F110W image. Its NIR 
photometric properties are very similar to those of the 



MYS O in 30 Dor (jRubio et all Il998t Brandner et al 
2001), indicating large luminosity. The source Y is lo- 
cated in a region of very large extinction, close (2") 
to a filled H II region (radio peak B) and the molec- 
ular cloud core. Thus, this area has all the necessary 
ingredients to be a star forming region. 



5 Summary 

In this contribution we have presented the first high- 
spatial resolution NIR images of NGC 604, obtained 
with the NIC2 camera aboard HST. These new NIC- 
MOS broadband images reveal new NIR point sources, 
clusters and diffuse structures. We found an excel- 
lent spatial correlation between 8.4 GHz radio contin- 
uum and 2.2 /im nebular emission. We also found that 
MYSO candidates appear aligned with those radio peak 
structures, reinforcing the idea that those areas are 
star-forming regions. Three different SOBAs can be 
recognized in the new NICMOS images, moreover, the 
brightest NIR sources in our images show properties 
that suggest the presence of RSG stars. 

This preliminary analysis of the NICMOS images 
shows the extreme complexity of the stellar content in 
the NGC 604 nebula. 
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